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Influence of disorder in the concentration of magnetic ions on the ferromagnetic 

phase transition in diluted (III,Mn)V semiconductors is investigated analyticaUy. The 
regime of smaU disorder is addressed, and the enhancement of the critical temperature 
by disorder is found both in the mean field approximation and from the analysis of 
the zero temperature spin stiffness. Due to disorder, the spin wave fluctuations around 
the ferromagnetically ordered state acquire a flnite mass. At large charge carrier band 
width, the spin wave mass squared becomes negative, signaling the breakdown of the 
ferromagnetic ground state. 



I. INTRODUCTION 



The prospects of application in spin electronic devices inspired much experimentaH and theoreticalf 
^ activity in studying the itinerant ferromagnetism in (HI,Mn)V diluted magnetic semiconductors. The 
C^eoretical efforts led to much understanding of the basic features of that material system, providing a 
C^ean-field descriptionj-and reasonable estimates of the ferromagnetic phase transition temperature from 
Itie spin- wave stiffnesaau . 

■ O There are two important issues, however, that have not been completely understood so far, namely, the 
!)le of disorder and interactions between the mobile charge carriers (holes) in the ferromagnetic phase 
ransition. In this paper, we address the issue of disorder in the concentration of magnetic Mn-dopants. 

le role of disorder is two-fold. On one hand, the variation in the Mn-concentration creates regions with 
local concentration of Mn-spins higher than the average. These regions tend to be magnetically polarized 
itinerant carriers at higher temperatures than the phase transition temperature for the homogeneous 
, J^^stem. The latter should lead to higher transition temperatures in the disordered system. On the other 
^*and, the disorder leads to localization of the mobile charge carriers, thus hampering the onset of long 
^■^ynge magnetic correlations, which would lower the transition temperature in the disordered system. 
0©resumably, the first effect prevails, if the localization of holes is weak, whereas in the case of strongly 
O^cfalized carriers (strong disorder) the second effect takes over. The enhancement of the transition 
'"temperature int, disordered systems has been reported in recent theoretical investigations bjt,^. N. Bhatt 
^d M. BerciucI, and J. Schliemann et aB. The results of this paper qualitatively agree withtra. However, 
e region of strong localization of the charge carries, where the suppression of the critical temperature 
esumably takes place is not addressed here. 



II. THEORETICAL MODEL AND APPROXIMATIONS 



o 

The microscopic description of the itinerant ferromagnetism in (III,Mn)V semiconductors is given by 
^e Kondo lattice Hamiltonian. In contrast to the Kondo system, in the system under consideration 
*^^(ie concentration of the charge carriers (mobile holes in the valence band) p is much lower than the 
''jCbucentration of the Mn-spins Nj^in- Thejiattcr rules out the Kondo effect. Moreover, the pd exchange 
Coupling J is larger than the Fermi energjoH. As a result, in the ferromagnetic state, the charge carriers 
(holes) are completely spin polarized. 

A hole magnetically polarizes the Mn spins in its vicinity, thus forming a magnetic polaron. Close 
to the ferromagnetic phase transition there is a large correlation length in the Mn spin system, and the 
dynamics of mobile holes is much faster than the dynamics of Mn spins. In this case, the magnetic 
polarization cloud created by a hole stays for some relaxation time tmu after the hole has left the region. 
The linear size of the polaron cloud can be estimated as i? vtmu , where v is the velocity of the hole 
(the motion of the hole inside the polarized region is assumed to be ballistic) . The condition of low hole 
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concentration, p <^ iVMn jcrB implies that there is no more than one mobile hole polarizing the cloud at a 
given moment of time. 

Therefore, the physical picture close to the phase transition point is that of slowly relaxing Mn-clusters 
(or polaronic clouds) polarized by the spins of mobile holes. The hopping of holes establishes long range 
correlations between the clusters. Due to the randomness in the Mn concentration, the number of Mn- 
spins in each cluster is randomly distributed around some average value. This picture is conveyed in the 
mathematical model below. 

The model consists of mobile holes moving on a lattice with sites i. The holes are represented by the 
fermionic operators a^, dia- The other ingredient is a system of Mn-spins 5/2 distributed randomly ewer 
a lattice with sites j. Since the concentration of Mn spins is larger then the concentration of holesEJ-H, 
there are many Mn-sites within a unit cell of the hole lattice. A hole polarizes Mn-spins in its vicinity. 
Because of the randomness in the Mn distribution, the number of polarized spins around the site i is 
random. 

The Hamiltonian of the model writes 

H = -t + h.c.) + J^^s.Sj,. (1) 

<i,i'> i j — ^ 

The first term describes the nearest neighbor hopping of holes. In this papet, a simplified one-band 
hole dispersion is used instead of the six-band model relevant for (Ga,Mn)As.Ll Whereas the one-band 
model is insufficient to give a quantitative prediction for the critical tempeinature, it reproduces the basic 
qualitative features of the hole system, including the effects of disorderoQ. In the interaction term, 

~ '^ta'^apo.ifj describes the spin of a hole on site i, and the operators Sji describe the Mn-spins in 
the vicinity of site i. The interaction term is written under the assumption of a rectangular hole wave 
function, so that the constant J does not depend on the distance between the hole and the Mn-spin 
within the interaction range. 

There is a total of Ni Mn-spins polarized by the hole at site i. The number Ni is taken to be random 
gaussian distributed with mean Np and variation A/y/Np 



p.,M NpiN-Np)^ 
P(Ni) cx exp 



2A2 



(2) 



In the experimental system the average Np is usually smallB, Np ^ 3. Therefore, the Poissonian distri- 
bution is actually more appropriate to describe the disorder in N. We take the Gaussian distribution 
for the sake of technical convenience. The case of the Poissonian distribution of Ni is left for future 
investigations. To restrict the influence of the unphysical region of negative Ni, the condition A ^ Np is 



adopted in numerical evaluations below (see Section III) 



To deal with disorder, we employ the replica-trick. In the replicated partition function for the Hamil- 
tonian ([l]), the interaction term is represented by exp[— J^- J2^=i J2a=i lo ^'''^i'^^ji]' where a denotes 
the replica-index, and r denotes the imaginary time. We decouple the interaction term by the Hubbard- 
Stratonovich transformation, introducing decoupling fields for the magnetization of Mn-spins and hole- 
spins. After the decoupling, the interaction part of the partition function assumes the form 



exp 



rP v-QTYaT pP 

e/ '^-e^v^+e/ d^Y. {^i^^i-^ti-n 

" •'^ <ij> a "^0 <i,j> 



(3) 



The term V , , XfJSf^ describes the Mn-spins in the cluster around site i under the influence of the 
^ — ''J 

effective field created by mobile holes. is the effective field created by a hole at site i and acting on 
the Mn-spin on site j. For a rectangular form of the hole wave function, we replace XfJ ~ Xf^ Vj = 
l,Ni. In what follows we use the Holstein-Primakoff representation for the total spin of the clusterB, 
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(4) 
(5) 
(6) 



where 6^, hi are bosonic creation and annihilation operators, describing the excitations around the ordered 
state of the cluster. There is a constraint on the number of bosons in each cluster {bf bi) < 2NiS, S ~ 5/2. 
The further calculations are carried out in the quadratic approximation in the bosonic fields bf ,bi , so 



we replace in the square root y 2NiS — bj'bi ~ ^J2NiS. 

Finally, under the assumption of the rectangular form of the hole's wave function, we approximate the 
effective field exerted by the Mn-cluster on the hole at the site i, X^j^i Yf/s"'^ « A^iYf^sf"^. Under the 
approximations above, the replicated partition function of the model can be written as 



/ D[aZ\D[aZ]D[bZ]D[bZ]D[X.r\D[Yr]e^Y> 



EE 



dr ■ 



-aZ{-dr + t^)aZ 



+br{-dr-xtiw)+t J2 EE(«" 



Jo 



'^A V2SN:,{x^;br + x^br) - (F,7(ara-ar) + yr(ar^+ar)) 

. 1 Jo ^ 



(7) 



Here afj, afj denote the Grassman fields associated with the mobile holes, and bf^, bf^ denote the complex 
variables associated with the bosonic operators bf,bi. We introduced the spinors a°'^,a°'^, which are 
formed as follows, af^ = «f ,a^p^, = (af^)+. 

The partition function (|^) is analyzed on the mean field level. The mean field approximations include 
neglecting the transverse fluctuations of spins (only z-components of fields are considered), and a static 
and spatially homogeneous ansatz for the fields Xf/ « X", Fj""^ « The averaging over the disorder 
with the distribution (^) results in a term in the effective free energy that mixes different replicas and 
looks as 



2Nr 



n 



(8) 

I a, 6=1 

The structure of term (|^) is typical for the replica treatment of the Anderson localization problem!. 
Under the mean field approximations made above, the strength of the disorder in the effective localiza- 
tion problem is proportional to the mean field value of the squared magnetization in the Mn-clusters 
Y°-Y^ . Therefore, expression (||) accounts for the localization of holes by the random potential from the 
frozen Mn-dusters. The treatment of repeats the nonlinear sigma- model approach for the Anderson 
localization^. The term is decoupled in the following fashion 



exp 



'2N, 



1 " r(3 

-EE / rfrdr'(ara^ar)(ara^ar') 



aCJ^ a(4z cxp 



i a, 6—1 k 



(9) 
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where Q"^'**^ is a matrix in the spin space with elements (Q"^'''^ )itit'- Here we concentrate only on the 
mean field solution for the matrix Q. After Fourier transformation with respect to imaginary times r, t', 
the mean field ansatz reads 

Qr-'"' = Ksign{uj)5^^>6ab. (10) 

Neglecting the transverse spin fluctuations, we assume a diagonal structure of the matrix A^^ in the spin 
sector 

Finally, the mean field partition function acquires the form 



a—1 UJn 



.A' 



2 



6=1 

2 



p 



6=1 



n N NN ^ 

a—l iOn i—1 a — 1 u 

+l3NNp^{SX'' -X''Y''/J} + -t^ ^ {S'2A"A'' + A'^A^'r^yVj^ 



a, 6=1 



-2SX''X''Y''IJ]] . (12) 

Here the sum over the frequencies aj„ denote the sum over the fermionic Matsubara frequencies for the 
fields a, a, and the sum over the bosonic Matsubara frequencies for the fields 6, h. N denotes the total 
number of sites i in the lattice. The exponent in ([l^ ) is quadratic in the bosonic (b, b) and the fermionic 
(a, a) variables, hence the bosons and fermions can be intergrated out. When integrating out the bosons, 
the restriction on the bosonic occupation number n;, < 2NpS is explicitely taken into account. 



III. MEAN FIELD EQUATIONS AND THEIR SOLUTIONS 

We adopt the replica symmetric mean field ansatz for the fields y° — Y, A° — X, A^'^ — AJ^^ Va = 1, n. 
In the replica symmetic approximation, the terms in (|l2|) involving the sum over different replicas give 
no contribution to the mean field equations in the replica limit n 0. 

Further, we change from the integral over wave vectors to an integral over the energies J d'^k... = 
J pd(f)dt, where Pd(() denotes the d-dimensional density of states (DoS). Under the above approximations, 
the following mean field equations have been obtained 

y_ J Y.IT'''^' (jA^P - Binh [/? {In, - X] 

Eli^''-^^/' cosh [/3 (1^,-0 A] • ^ ' 

Here Y is the average magnetization of a single Mn-spin in units of J. Mn-spins are magnetically polarized 
by the effective field of the holes X, which is reflected by the RHS of Eq. ([l3|). The average magnetization 
of a hole X (in units of J) is given by the equation 
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+JTy p{e)de{ — 



+ ^ - e + NpY + iY'^A'^^sign{uJn) 



iujn + /i — e — NpY + lY'^ K^^sign[ijJn) 



(14) 



Eq. (14) describes the magnetic polarization of itinerant carriers (holes) in the effective field NpY created 
by Mn-spins. Disorder enters Eq. ( pH ) through the disorder strength A and the fields A^^ which, on 
the mean field level, are proportional to the inverse mean free time for holes. At zero temperature, the 
inverse mean free time for holes with a spin projection a is given by I/to- = 2y^A^g.. Since only the 
scattering by frozen magnetic moments is taken into account on the mean field level, the influence of 
disorder disappears completely in the paramagnetic phase. In the paramagnetic phase, the scattering by 
dynamical magnetic fluctuations in Mn spin clusters affects the hole motion. The account for that effect 
is beyond the static mean field approximation adopted here. The analysis of that and other fluctuation 
effects beyond the mean field approximation is left for future investigations. 



Equations for the mean field values of A^^. can be written as 



A'; 



A2 

N,. 



w 



p{e)de- 



i sign{LUn) 



w 



^l-eT NpY + iY^A'^^ sign{uJn) ' 



(15) 



The upper sign in the denominator corresponds to a =| and the lower sign corresponds to a —I. For 
d = 2, 3, 1 replace the density of states /3d(e) in jisl ) by a constant average density of states pd{e) = p. 
The example of a spherically symmetric band e(k) = — M^cos(|k|) shows that for large W, the contribution 
of the singularity of the DoS near the band edge e = W is small in both two and three dimensions, hence 
the above replacement is justified. Upon that replacement, the integrals over the energy e in Eqs. (p^, 
can be performed explicitely, and the equations for X and A^-a assume the form 



X = JpTj2 log 

LJ„>0 



[iW-p + NpYf + K + Y^A^^^f][{W + P + NpYf + K + Y^^XO^] 



J 



F 



[{W + p- NpYf + K + y2A- )2][(^^ _ ^ _ NpYf + K + y^A- )2 
l' + (Au)'}- 



(16) 



A' 



TT 



arccos 



NpY - W - 



\^ \NpY -W~p + i{tOn + i^2Af^)| 
NpY + W-P 



\NpY + W-P + i{ujn + Y^A^^)\ 



UJn > 0, 



(17) 



A' 



N„ 



NpY -^W + p 



arccos 



\^ \NpY ~W + p + t{uj,, + y2Aj'^ 
NpY + W + P 



\NpY + W + P + z(c^„ + r2A- )| 



iOn > 0. 



(18) 



We analyze the mean field phase diagram by solving numerically mean field equations (|l|), (|l|)-([8|), 
taking the parameters close to those of (Ga,Mn)As. The exchange constant for (Ga,Mn)As equals Jq « 55 
meV nm'^. The constant Jq is adopted in the model with a (5-like interaction between a hole and Mn-spins 
of the fornia 



Hi, 



./o5IS/s,,5(r,-R/). 



(19) 



Comparison of (|9|) with the interaction term in (|l|) for the homogeneous case Ni = Np gives the relation 
J = Jo/Np. The number Np can be evaluated as Np ~ Nmu/p- For Nmu = 1-0 nm~'^, p = 0.35 nm~'^ 
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we infer Np ~ {Nmu/p) ~ 3. For the simplified band structure adopted here, the average density of 
states p and the width of the band W cannot be directly obtained from experimental data. The chemical 
potential is fixed by the condition of constant concentration of holes p, or constant filling factor v of the 
charge carrier lattice i of the model (El) . 

The mean field equations (|l|), - (|lll) were solved numerically for J = 10 meV, and different 



values of Np, p, W, and A. In all cases we found that the critical temperature increases with disorder 
strength A. The temperature dependence of the magnetizations of holes X and Mn-spins Y for the set 
of parameters J = 10 meV nm'^, Np = 5, W = 100 meV, p — 0.01 meV~^ nm~^, v = 0.25 and different 
values of disorder strength A is shown in Fig. |^. 

IV. SPIN WAVES IN A DISORDERED FERROMAGNET 

There are massless Goldstone modes in the ordered phase of a clean system - the spin waves. The zero 
temperature spin stiffness has been used in Refu to give an estimate for the critical temperature, which 
is much closer to experimental data than the mean field result in the regime of small concentration of 
charge carriers and large exchange coupling J (as compared to the Fermi energy). 

The spin wave excitations are described by bosonic fields b, b. In this section, we calculate the zero 
temperature dispersion of the spin-wave excitations around the saddle point ( |l3| ) - ( p^ described above, 
and use the result to analyze qualitatively the influence of disorder on the critical temperature calculated 
from the spin wave stiffness. Since the role of disorder increases at low temperatures, the calculation of 
the zero temperature spin stiffness using the mean field description of the disorder effects (the fluctuations 
of the field Q are not taken into account) is perturbative in the disorder strength A. It can be used only 
to describe the qualitative behavior of Tc. The zero temperature spin wave dispersion is given by 

n^Dq^+ml^. (20) 

The spin waves acquire a mass msw- Both D and A also have imaginary parts that describe the 
attenuation of spin waves in a disordered system. At zero temperature, the approximate expression 
for the spin stiffness reads 



3NpS V N^p. 

where we neglected the imaginary part and assumed np <C NpJ, which is consistent with the regime of 
low concentration of holes, vp denotes the Fermi velocity. One can see from expression ( |2l] ) that the 
spin stiffness grows with disorder strength A at small disorder. This in turn results iii-an increase of the 
critical temperature. The critical temperature can be estimated according to formulaci 

T, « Dkl (22) 

with the Debye wave vector kn = {Gn^Np)^/^ . 

The critical temperature versus the disorder strength A calculated in the mean field approximation is 
shown in Fig. |[ The critical temperature grows strongly with disorder. In contrast, the evaluation of the 
critical temperature from the spin wave stiffness (^) gives a much weaker dependence on the disorder. 

The zero temperature mass of the spin wave excitations msw is given by 

-Iw « ^^^-P(A?T - a;,) -^Jd- {(An)^ + (An)'} . (23) 

where we neglected a small imaginary part. The dependence of the spin wave mass 'rn'^^r on the hole 
bandwidth W at two different sets of other parameters is shown in Fig. ^. At large bandwidth W, rngy^ 
becomes negative (see Fig. ||), which signals the breakdown of the coUinear ferromagnetic order. The 
transition to negative rnj^y^ occurs at lower values of W for a smaller average number of spins Np in a 
magnetic cluster. 

Presumably, in the regime < the ground state is noncoUinear. The existence of a noncoUinear 

ground state in disordered (III,Mn)V semiconductors has been suggested recentlyQ. 
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V. SUMMARY AND DISCUSSION 



In this paper we investigated the influence of disorder on the critical temperature of an itinerant 
ferromagnet. The theoretical modeling in the regime close to the phase transition was based on the 
picture of magnetic polaronic clouds (magnetic clusters) with large relaxation time. The clusters are 
polarized by mobile holes hopping between them. The regime of small concentration of mobile holes 
considered here is relevant for the ferromagnetism in (III,Mn)V compounds. The disorder in the Mn 
concentration is naturally included as a random number of spins in each magnetic cluster. At the critical 
point, the clusters begin to percolate, developing infinite range magnetic correlations. 

We found that the critical temperature grows with disorder, both in the mean field approximation 
and if calculated from the zero temperature spin stiffness. However, the influence of disorder turns out 
to be much stronger in the mean field evaluation of Tc than by the estimation from the spin stiffness. 
In both approaches, disorder is taken into account perturbatively. In the mean field calculation, the 
perturbative treatment of disorder can be justified, if the critical temperature is high (which is the case 
in the experiment), and the localization is suppressed. 

In contrast, the estimation of Tc from the zero temperature spin stiffness is much less reliable. At zero 
temperature, the higher order localization corrections become important, and perturbative account for 
disorder may result in the underestimation of the disorder effect on the spin stiffness. The latter, in turn, 
leads to the underestimaiton of the influence of disorder on Tc. At the same time, the calculations of 
the spin stiffness correctly reflect the general tendency that the spin stiffness, and therefore the critical 
temperature, grows with disorder. l-l 

The physical reason for the increase of Tc at weak disorder was proposed in Ref.B. In the disordered 
system, a mobile hole spends more time in the regions with higher Mn concentration, thus the effective 
magnetic coupling increases. However, at strong disorder, the critical temperature should eventually 
decrease because of localization of the charge carriers that mediate the magnetic correlations. 

In this work we neglected the localization corrections that should make the motion of the holes diffusive. 
The work on this subject is in progress. 

We thank A. H. MacDonald, J. Schliemann, J. Konig, and K. Scharnberg for useful discussions. 
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FIG. 1. Relative magnetization of Mn-spins Y/{JS) (> 0) and holes —X/{2J) (< 0) versus temperature at 
different disorder strengths. The critical temperature increases with the disorder strength A. In order of increasing 
Tci A = 0.1; 1.0; 2.0; 3.0. Other parameters: J = 10 meV nm^, Np = 5, W = 100 meV, p = 0.01 meV^nm"^, 
V = 0.25. 
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FIG. 2. Relative change of the critical temperature with the strength of disorder rc(A)/Tc(A = 0) as calculated 
from the mean field equations (solid line), and from the zero-T spin wave stiffness (dots). J = 10 meV nm^, 
Np = 5,W = 100 meV, p = O.OlmeV^nm-^, z/ = 0.25. 
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FIG. 3. Zero temperature spin wave mass versus the bandwidth W. J = 10 meV nm^, Np = 5, ii = 0, p = 1.0 
meV"^ nm"^. Solid line: A = 0.5. Dashed line: A = 0.3. 
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